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Introduction 

This  paper  presents  an  algorithm  for  computing  steady-state  probabili¬ 
ties  and  mean  queue  lengths  and  delay  times  for  M/M/c  queues  with  two  priority 
classes.  Cobham  [2]  has  computed  mean  delay  times  for  multipriority  M/M/c 
queues  with  homogeneous  service  rates.  Heyman  [3]  pointed  out  the  need  to 
compute  mean  delay  times  where  service  rates  for  different  priority  classes 
are  unequal.  The  algorithm  presented  below  computes  expectations  for  a 
system  with  two  priority  classes  having  heterogeneous  rates;  in  addition, 
it  can  compute  steady  state  probabilities  for  individual  states.  However 
it  is  doubtful  whether  this  algorithm  could  be  modified  to  a  practical 
algorithm  for  systems  involving  more  than  two  priority  classes.  The 
preemptive  case  is  simpler  than  the  nonpreemptive  case;  an  algorithm  for 
it  can  be  developed  using  ideas  presented  in  this  paper. 

The  general  approach  to  this  analysis  is  a  computational  method  devel¬ 
oped  by  Neuts  [8],  The  simpler  problem  of  single-server  two  priority 
Markovian  queues  was  solved  by  Miller  [5].  The  approach  consists  of  parti¬ 
tioning  the  full  state  space  into  blocks  and  discovering  special  structure 
of  the  invariant  probability  vector  in  relation  to  these  blocks.  It  is 
assumed  that  the  reader  is  familiar  with  [5, 6, 7, 8]. 

In  addition  to  solving  the  problem  posed  by  Heyman,  this  algorithmic 
approach  is  of  interest  because  it  explores  a  new  and  complicated  application 
of  several  ideas  from  the  field  of  computational  probability:  quasi  birth 
and  death  processes,  M/G/l  paradigms,  and  matrix-geometric  invariant  vectors. 

Using  single  precision  arithmetic  the  algorithm  appears  to  work  well 
for  low  and  moderate  utilization  factors,  but  there  is  some  deterioration 
in  the  calculations  for  high  utilization.  This  application  demonstrates 
the  need  for  developing  good  error  analyses  for  this  type  of  computation. 

This  application  may  be  a  good  vehicle  for  developing  some  such  analyses. 
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State  Space  and  Transition  Matrices 

Consider  an  M/M/c  two  priority  nonpreemptive  queueing  system  with 
arrival  rates  A^  and  am*  service  rates  and  •  The  state 

space  can  be  described  as  follows.  Let  x.  .  ,  be  the  state  with  i 

i » .1 » 

first  priority  customers  waiting,  j  second  priority  customers  waiting, 
k  first  priority  customers  in  service  and  c-k  second  priority  customers 
in  service,  i,  j>0,0<k<c.  From  them,  define  the  blocks 


H  =  (x  ,  jO  £  k  £  c}  ,  i,j  >  0  . 

Let  x  be  the  state  with  m  customers  in  service  of  which  n  are 
m,n 

first  priority  customers  and  no  customers  are  waiting  in  queue, 

0  n  <  m  <  c  .  Then  define  the  blocks 

H  =  {x  |0<n<m},  0<m  <  c  . 

m  m,n'  —  —  —  — 

(Note  that  H  =  H_  .  .  )  Thus,  the  state  space  is 
c  U,  u 

S  .  (  U  H  \  U  (  U  U  H  \  . 

\  m=o  m/  \i^>o 


Now  consider  the  probabilistic  transition  rates  for  the  non¬ 
preemptive  M/M/c  two  priority  system  on  this  state  space.  The 
matrices  of  transition  rates  between  the  blocks  defined  above  can  be 
denoted  as 


H 

H 

H 

H 

H 

H 


i— Hi+1  • 
i— Hi-1  * 

iTrvi.j* 

iTT^Hi,j+l’ 

i7T*Hi-i,r 

oTT’Vj-r 


0  <  i  _<  c  — 1 
1  <  i  <_  c 
i,j>  0 
i,j-0 

ill,  jiO 
j  1  1  • 


-  2  - 
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All  transitions  correspond  to  arrivals  and  departures  thus  there  are  no 
transitions  within  individual  blocks.  Therefore  the  submatrices  of  the 
transition  rate  matrix  corresponding  to  transitions  within  a  block  are 
all  diagonal: 


-D.  : 

i 

H.  - » 

l 

H±  , 

0  < 

i 

<  c 

-D  : 
c 

H.  t— * 
l.J 

Hi,j  * 

i,  j 

> 

0 

Ihe  transition  structure  is  depicted  in  Figure  1. 

The  matrix  K.  is  (i+1)  *  (i+2)  dimensional;  (K.) .  equals 
i  i  j  t  k 

A^  for  k  *  j+1  ,  for  k  =  j  ,  and  0  otherwise  (1  <  j  <  i+1  , 

1  <  k  <  i+2)  .  The  matrix  J.  is  (i+1)  x  i  dimensional;  (J.) 

~  ‘  1  i],lt 


equals 

(i-j+l)y2  for  k  =  j  , 

(j-l)y1  for 

k  = 

j-l  ,  and  0 

otherwise 

‘1™* 

V  1 

<  i+1  ,  1  <  k  <  i)  . 

The  matrices 

L1  ’ 

L2  ,  ,  and 

M2  are 

(c+1)  x 

(c+1)  dimensional; 

equals  A^I  ; 

equals  A2I  ; 

equals  (c-j+l)y2  for  k  *  J+1  ,  (j-l)y  for  k  *  j  ,  and  0  otherwise 
(1  <  j  <  c+1  ,  1  <  k  <  c+1)  ;  k  eclua^s  Cc— j-+l) yi2  for  k  =  j  , 


(j-l)y^  for  k  *  j-1  ,  and  0  otherwise  (1  <  j  <  c+1  ,  1  <  k  <  c+1)  . 

The  matrix  D.  is  (i+1)  x  (i+1)  dimensional;  (D.).  equals 
1  1  J  »k 

(j-l)y^  +  (i-j+l)y2  +  A  for  k  =  j  and  0  otherwise  (1  <  j  <  i+1  , 

1  <  k  <  i+1)  . 

The  blocks  H  ,  0  <  m  <  c-1  ,  and  H.  .  ,  i,j  >  0  ,  can  be 

m  --  - 

combined  into  super  blocks: 

!_1  =  CJ  H  .  !.=  U  H  ,  i  >  0  . 

i=o  j=o 

-  3  - 
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The  states  in  1_^  are  exactly  those  corresponding  to  the  existence  of 

Idle  servers.  The  states  in  I  ,  i  >  0  ,  are  those  wfth  exactlv  1 

i 

first  priority  customers  awaiting  service.  The  matrices  of  transition 
rates  between  and  within  these  blocks  can  be  denoted  as 


B-i,  -  r 

^l"— 

x-i 

T-l— 

To 

Bo,-i: 

*0  — 

T-1 

ttf 

O 

O 

T  ^ 

To 

^0 

o 

< 

Xi+1  ’ 

i  > 

A  * 

I  » 

i  > 

V 

i 

Ji  ’ 

A2: 

1  i 

Ti~i  * 

i  > 

These  transitions  are  depicted  in  Figure  2.  The  above  matrices  of 

transition  rates  can  each  be  partitioned  into  submatrices  corresponding 

to  transitions  between  the  subblocks  of  I.  ,  i  >  -  1  .  These  sub- 

i 

matrices  were  defined  earlier.  The  number  of  submatrices  in  one  of  these 

matrices  varies  according  to  the  number  of  blocks  in  the  corresponding 

superblocks.  (B  ,  .).  ,  equals  K.  .  for  k  =  j+1  ,  -  D.  for 

-1,-1  J,k  J-l  j-1 

k  =  j  ,  J  ^  f°r  k  *=  j-l  ,  and  0  otherwise  (l<j<c,l<k  <  c)  . 

(B  -  .),  ,  equals  K  for  (j,k)  =  (l,c)  and  0  otherwise 
J > K  C“X 

(1  <  j  <°°,  1  <  k  <  c)  .  <Bn  ).  ,  equals  J  for  (j,k)  =  (c,l)  and 
0  otherwise  (1  <  j  <  c  ,  1  <  k  <  °°)  .  (B^  ^  equals  L ^  for 

k  =  j+1  ,  -  Dc  for  k  =  j  ,  for  k  =  j-l  ,  and  0  otherwise 

(l<j<w,l<k<°°)  .  (An)  .  equals  L.  for  k  =  j  and  0 


-  5  - 
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Partitioned  state  space  (superblocks)  and 
transition  scheme. 


_  6  _ 
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otherwise  (1  <  j  <  <»,  1  <  k  <  °°)  .  (Ad.  .  equals  L  for  k  =  j+1  , 

“  "  ±  J  9  - 

-  D  for  k  =  j  ,  and  0  otherwise  (l<j  <aj,  1  <  k  <  °°)  .  (A„)  .  , 
c  -  -  ^  J  >  *■ 

equals  for  k  =  j  and  0  otherwise  (1  <  j  <  °°,  1  <  k  <  °°)  . 


Using  this  block  structure  and  notation,  the  transition  rate  matrix  for 
the  M/M/c  two-priority  nonpreemptive  queueing  system  is 


The  process  described  above  is  a  quasi  birth  and  death  process 


L 4 , 8  1 .  Consequently  its  steady-state  probability  vector  is  of  matrix- 
geometric  form  [4,8] : 


(  2-1 


2o 


V 


2oR 
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The  Matrix-Geometric  Rate  Matrix 

The  rate  matrix,  R  ,  of  the  quasi  birth  and  death  process  is  the 
minimal  solution  of 

AQ  +  RAj  +  R2A2  =  0  .  (2) 

Furthermore,  from  the  block  structure  of  the  process  and  the  interpretation 
of  the  rate  matrix  [5,7,8],  R  must  have  the  structure 


where  each  submatrix  has  dimension  (c+1)  x  (c+1).  Substituting  (3)  into 
(2)  gives  a  system  of  equations: 


(A) 

(5) 


8 
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These  can  be  solved  numerically  as  follows.  First  consider  (4).  Let 


L  =  D  LL.  ,  M  =  D-1M.  ,  R.  =  D"V.D  ; 

l  cl  1  cl  0  cOc 


then  (4)  becomes 


'V,  'V  O/  'V 

R0  =  L1  +  Ro  M1 


which  can  be  solved  using  the  usual  iterative  approach  useful  for  matrix- 
geometric  rate  matrices  [8]:  Let 

%  a, 

so  =  0  *  si  =  Li  +  so  Mi  *  »  si+i  =  Li  +  si  Mi  • 

Then  J  ,  termwise.  The  usual  procedure  is  to  continue  iterating 

until  the  maximum  termwise  difference  between  successive  iterates  is 
smaller  than  some  £  ,  e.g.  £  =  10  ^  . 

Now  consider  (5);  it  is  equivalent  to 

'Xj  *Xj  ^  i—  l'V'X/  'V  %  s\j  'Xj  'V  'Xj 

ri  •  R1-1L2  +  jSl  RjRi-,Ml  +  <R0Ri  +  RiV  M1  •  1  i  1  (7) 


For  i  =  1  ,  R^  can  be  found  by  using  the  above  solution  for  and 

then  using  a  similar  iterative  procedure  which  will  converge  monotonically 
to  R^  .  For  i  >  1,  continue  recursively,  using  the  solutions  for 

R.  ,  j  <  i  ,  from  previous  steps  and  using  the  iterative  procedure  to 

n, 

get  R.  .  The  desired  R.'s  are 
&  l  l 


R.  =  D  R.D  ,  i  >  0  . 
1  c  l  c 


Thus  the  rate  matrix  R  can  be  computed  up  to  any  level  of  truncation. 
In  this  study  the  computation  of  the  R^'s  was  truncated  when  the 

iterative  procedure  stopped  in  the  first  iteration,  in  which  case  the 
value  was  set  to  0  . 
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There  is  an  internal  accuracy  check  which  can  be  used  in  the  above 
numerical  computation.  Let 

a, .  0°  'K, 

R  =  Z  R  . 

i=0  1 

From  (6)  and  (7)  it  follows  that 

%  %  l\,  <\j  2 

R*  =  L  +  R*  L2  +  (R*)  M  .  (8) 

'V* 

The  usual  iterative  method  can  be  used  to  solve  (8)  for  R  .  This  matrix 

n, 

can  then  be  compared  to  the  sum  of  the  individual  solutions  R^  ,  i  >  0  . 
This  check  was  performed  and  virtually  no  error  detected. 

The  M/G/l  Paradigm 

According  to  the  theory  of  matrix- geometric  invariant  vectors,  if 


0  =  TrPg  (9) 


then 


7t  —  (  »  3Iq  »  *  33 2  *  •  •  • ) 


where 

and 


corresponds  to  the  vector  on  the  super  block 


I. 

l 


i  > 


Ii+1 


=  TT  .R 

-X 


i  >  0 


This  can  be  used  with  (9)  to  temporarily  reduce  the  problem  of  solving 
(9)  to  consideration  of 


(0  0)  = 


(  33.!  330) 


-  10 


T-440 


whore 


is  a  matrix  with  negative  diagonal  entries,  nonnegative  off-diagonal 
entries,  and  row  sums  equal  to  0  ;  thus  it  can  be  thought  of  as  a 
transition  rate  matrix  for  a  Markov  process  on  I  •  Letting 


c0  ■  M2  '  Cl-  -°c  +  R0M  •  C2 


L-  +■  RM,  ,  C  -  R.  M  ,  i  >  3  , 
2  11  i  l-l  1 


gives 


Note  that  this  has  the  M/0/1  paradigm  structure  discussed  by  Lucantoni 
and  Neuts  [4].  Following  their  approach,  it  will  be  easier  to  work  with 
an  embedded  Markov  chain.  Define  the  submatrices  of  transition  probabil 
ties  for  this  chain  as 

%  -1  ^  -1 

K,  =  D.aK.  ,  J.,  =  D/J.  ,  0  <  i  <  c  , 

i  ii  i  ii  -  - 

^  _2 

C,  =  D  C,  +  I  , 

1  cl 

C  =  D~V  ,  i  =  0,  2,  3,  ... 
i  c  i 
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This  is  a  Markov  chain  on 


I  .  UI 

-l 


0 


c-1 

u 

m=0 


m 


U 


U 

j=0 


0,j 


Each  transition  of  this  chain  corresponds  to  a  transition  of  the  original 

process.  Each  transition  of  the  original  process  out  of  a  state  in 

\ 

I  ^  UIq  a^S0  a  transition  for  this  chain.  If  the  process  made  a 
transition  from  1^  to  1^  ,  then  the  chain  will  make  a  transition 
from  Iq  to  Iq  ,  the  target  state  being  the  first  return  state  in 
1^  .  Thus  the  chain  can  make  transitions  from  states  of  1^  into  them¬ 
selves. 

The  following  analysis  uses  ideas  from  Lucantoni  and  Neuts  [4]. 

n, 

Consider  the  Markov  chain  PT  _  on  this  state  space.  Let  0  be 

'-l^o 

the  (c+1)  x  (c+1)  matrix  of  hitting  probabilities  of  states  of  H  . 

» 1 

starting  at  H  . ,  j  >  0  .  G  must  satisfy 
U  >  J"’"! 

oo 

G  =  Z  C  GV  , 

V=0 

and  can  be  found  iteratively  starting  with  an  initial  0  matrix;  see 

[4].  Let  V  be  the  first  passage  probabilities  of  hitting  states  in 

H.  starting  from  states  in  H-  .  with  .  a  taboo  set;  then 

O.j+l  0,j  0,j 


V 


00  f\j 

v=2 


„v-2 


Let  W  be  the  first  passage  probabilities  of  hitting  states  in 

with  H.  .  a  taboo  set  starting  from  states  in  IL  ,,,  ;  then 

0,j  0,j+l 


H 


0,j+l 


W 


z 

v=l 


'V 

Cv  G 


v-1 
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a. 

Finally,  let  be 

visits  to  states  in 

H  .  ;  then 
9  J 


the  matrix  whose  entries  are 

...  from  states  in  . 

0,j+l  0,j 


the  expected  number  of 
before  returning  to 


a, 


=  I 
v=0 


Using  the  above  relationships 


vwv  =  V(I-W)  1 

n, 

can  be  computed. 


Probabilities  of  Idle  Servers 


Now  consider  an  invariant  vector  for  the  Markov  chain  with  state 

n, 

space  I  UIn  and  the  Markov  transition  matrix  P  : 

1  U  Q 


(5-l  ’  V  =  (5-i  *  20>  Pl_1ui0 


(10) 


Partition  the  vector  over  the  blocks  H  ,  0  <  m  <  c-1  ,  and 

m  -  - 

Vi  *  1  - 0  : 


(--i  »  v  =  (2-i,0  ’  --1,1  ’ 


’  --l,c-l  *  -0,0  ’  -0,1  ’  -0,2  ’  m-m) 


From  the  structure  of  this  process  and  an  important  property  of  taboo 
probabilities  (Theorem  1  of  [5],  or  see  Chung  [1],  p.  53),  it  follows 
that 

-0,1  =  -0,0  S1 


This  result  applied  to  (10)  gives 
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The  dimension  of  the  above  square  matrix  is  (c+1)  (c+ 2)/2.  For  a 
moderate  number  of  servers  c  ,  the  invariant  vector  (z  ,  z_  ) 

I  k)  ,  vj 

can  be  computed  using  existing  numerical  techniques.  The  corresponding 
invariant  vector  for  the  process  has  component  vectors 


y  .  .  •-  2  D.1,  0  <  i  <  c-1 

—  l,i  —l,i  1  -  - 


^0,0  -0,0  Dc 


-1 


In  order  to  get  a  normalizing  constant  to  convert  this  into  the  invariant 
probabilities,  recall  that  the  proportion  of  idle  servers  must  equal 


P  =  1  - 


\  +  S 


CP, 


CP, 


where  p  is  the  utilization  factor.  Thus,  let 

c-1 


0.-1 


i=0 


Z-l.i  • 


Then  the  invariant  probabilities  on  I  ,UH  can  be  computed: 

—  1  U  *  J 


2  ,  ,  =  y  ,  .  (l-P )/€  ,  0  <  i  <  c-i  , 


-o,o  ~  ^0,0  • 


Probabilities  of  States  with  Customers  Waiting 

Now  it  is  possible  to  build  the  state  space  back  up,  computing  the 
invariant  probabilities  of  additional  states.  In  order  to  compute  the 
invariant  probability  vector  (tt^  ^  »  Iq  2  ’  ^0  3  *  *  * "  ^  over 


OO 
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it  is  necessary  to  depart  from  the  approach  of  Lucantoni  and  Neuts  |4) 

n.  % 

because  is  singular.  Instead,  consider  S.  ,  the  matrix  whose 

entries  are  the  expected  number  of  visits  to  states  in  ...  from 

0,.]+i 

states  in  H_  .  before  hitting  ,  ,  k  <  j+i  ; 

0,j  b  0,k 

%  oo  n.  v_i_i 

S.  =  E  C  G  (I-W) 

1  V=i+1  V 

This  follows  by  a  simple  sample  path  argument  similar  to  the  derivation 

n, 

of  Sj  earlier.  The  invariant  vector 


(5-l  ’  50,o 


-0,1  ’  -0,2 


...) 


for  the  Markov  chain 


a. 


UI 


0 


must  satisfy 


-0,1+1 


l 

=  E 

j-0 


n, 

zn  .  S.,_  .  ,  i  >  0 

-0,3  i+l-j 


This  is  a  special  case  of  the  fundamental  result  for  taboo  probabilities 
(Theorem  1  of  [5]  or  see  Chung  [1],  p.53).  The  invariant  probabilities 
for  the  process  must  therefore  satisfy 


...  =  E  TT_  .  D  S.,,  .D 

~°’1+1  j=0  -°.-l  C  1+1-3  C 


-1 


i  >  0  . 


These  can  be  computed  recursively,  starting  from  q  which  has  already 
been  computed. 

Finally  it.  .  ,  i  >  1  ,  j>0,  can  be  computed  using  the 

-1 » J  ~  ~ 

matrix-geometric  structure: 


-i+1.3  k^0  j~k 
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Thus  one  can  compute  the  invariant  probability  vector  to  any  level  ot 
truncal  ion  0  <"  i  <  I  ,  0  <  i  <  J  . 

Sums  and  Means 

By  summing  the  above  probabilities  it  is  possible  to  get  a  separate 
calculation  of  the  probability  of  no  idle  servers  (this  can  be  used  as  a 
consistency  check  on  the  numerical  calculation)  and  the  mean  number  of 
each  type  of  customer  awaiting  service  (then  Little’s  formula  can  be 
applied  to  compute  mean  delay  for  each  class). 


£  it.  ,et  =  n  (I-R*)  1  et 
i»  j-0  -l’J“  - 


(ID 


00  °°  «-  *  *  -2  t 

E  i  Z  IT.  .<3-  =  n  R  (I-R*) 

i=l  j=0  -1,1_ 


(12) 


oo 


-1  it  re  ^  ^  Ct_c  W)  m 


E  j  E  it .  =  A"x  II  (Svt'(I-Sw/)  -  S"")  M  « 

j=l  i=0  1,J“ 


(13) 


where  e*"  is  a  (c+1)  dimensional  column  vector  consisting  of  all  I's  , 
and 


n 


10.0  <-(0))'1 


(0)  00  ^  -1 

>  '  =  D  E  S,  D  i 

C  k=l  k  C 


.(1) 


oo  ^1 

D  E  k  S,  T) 

C  k=l  k  C 


Computational  Experience 

The  algorithm  described  above  has  been  programmed  in  single-precision 
Fortran  and  run  on  GWU's  IBM370/3031.  Cases  with  c=5  were  run.  Seventy 
five  different  cases  were  run  corresponding  to  all  combinations  ol  p  -  .2  , 
.  J  ,  . 8  ,  A1/a2  -  .23  ,  .5  ,  1,2,4,  and  =  -25  ,  .5  ,  1,2,4.  Executio, 

times  for  computing  state  probabilities  and  expectations  varied  from  approx¬ 
imate  lv  10  seconds  per  case  with  P  “  .2  to  approximately  60  seconds  per  case 
with  .  ■  .»  . 
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Various  consistency  checks  were  used  in  the  computational  prucoUur. 

t)  independent  calculations  of  R  and  Z  R.  were  compared;  ii)  t  u 

i=0  1 

row  sums  of  C.  were  compared  with  unity;  iii)  the  row  sums  of 

't  a.  n.  n, 

+  C  +  S|  Cq  were  compared  with  unity;  iv)  the  total  probability 

computed  was  compared  to  unity;  and  v)  moments  were  computed  direct Iv 
from  the  state  probabilities  and  compared  with  values  computed  from 
equations  (11)  and  (12).  These  consistency  checks  generally  agreed  to 
5  or  more  digits. 

The  mean  delay  for  cases  with  homogeneous  service  rates  were 
calculated  using  Cobham's  [2]  approach.  The  values  agreed  with  those 
computed  by  the  above  algorithm  except  in  the  case  p  =  .8  where  a 
discrepancy  appeared  in  the  fourth  digit.  This  case  illustrates  the 
need  for  a  more  complete  error  analysis  in  this  type  of  calculation. 
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Relative  Mean  Delays 
(5  servers,  .2  utilization) 


Xl/X2 


V^2 


ill 

.50 

1.0 

2.0 

.25 

.0003698 

.001947 

.0004961 

.002613 

.0006405 

.003333 

.00007681 

.003922 

.50 

.0004947 

.001251 

.0005728 

.001453 

.0006720 

.001705 

.0007704 

.001948 

J_.J) 

.0007980 

.0009978 

.0008211 

.001026 

.0008514 

.001064 

.0008841 

.001105 

2.0 

.001454 

.0009104 

.001394 

.0008738 

.001308 

.0008214 

.001209 

.0007597 

4.J) 

.002846 

.0008955 

.002665 

.0008431 

.002384 

.0007599 

.002020 

.0006490 

Table 

1.  Values 

of  and 

P2W2  for 

c— 5  9  P  ”  •  - 

and  X 

j/X2  =  .25  ,  . 

5  ,  1,  2,  4 

,  and  n  /i!0 

i  *- 

,.5,1,  2,  4. 


1  <3 


4.0 


.0008479 

.004441 


.0008')  10 
.0021 14 


.000912  1 
.  00 1  1 4 1 


.001118 

.000702! 


.001650 
.u005  516 


.25 
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Relative  Mean  Delays 
(5  servers,  .5  utilization) 


.25 

VX2 

.50 

1.0 

2_.  0 

4.0 

.01488 

.02195 

.03082 

.03892 

.04414 

.25 

.1474 

.2136 

.2854 

.  3425 

.  3782 

.01844 

.02286 

.02922 

.03636 

.04250 

L 

to 

.07343 

.09859 

.1258 

.1540 

.  1765 

.02898 

.03128 

.03476 

.03911 

.04  146 

1.0 

.05795 

.06257 

.06953 

.07822 

.08690 

.05220 

.05224 

.05226 

.05226 

.05223 

2  J3 

.05202 

.05192 

.05172 

.05168 

.05165 

.1011 

.09835 

.09364 

.08669 

.07809 

4.0 

.05083 

.04964 

.04759 

.04442 

.04030 

Table  2.  Values  of 

ui»i 

and  iJ2W2  for  c=5  *  P  =  ‘  ^  ’ 

VX2  = 

.25  , 

.5  ,  1,  2,  4,  and  Pj/l^  = 

.25  .  .5  ,  1,  2,  4  . 
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Relative  Mean  Delays 
(5  servers,  .8  utilization) 


xx/x 


2 


.25 

.50 

1.0 

2.0 

_._25 

.2113 

3.9591 

_J0 

.1159 

1.3541 

.1726 

1.8911 

.2437 

2.6118 

1.0 

.1319 

.6596 

.1511 
.  7553 

.1846 

.9218 

.2405 

1.1956 

2.0 

.2462 

.5983 

.2663 

.6369 

.2976 

.6880 

4.0 

.4589 

.5520 

Table  3.  Values  of  and  ^2W2 

for  c=5  , 

p  =  .5 

A1/A2  =  .25  ,  .5  ,  1,2,  4, 

and  y  /u2 

=  .25 

.5 


2,  4. 


.  0 


.  3052 
1,5113 
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To  cope  with  the  expanding  technology,  our  society  must 
be  assured  of  a  continuing  supply  of  rigorously  trained 
and  educated  engineers.  The  School  of  Engineering  and 
Applied  Science  is  completely  committed  to  this  ob¬ 
jective. 


